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9 Abstract 

We investigate the dependence of the amount of the observed galactic cosmic 
ray (GCR) influx on the solar North-South asymmetry using the neutron 
count rates obtained from four stations and sunspot data in archives spanning 
six solar cycles from 1953 to 2008. We find that the observed GCR influxes at 
Moscow, Kiel, Climax and Huancayo stations are more suppressed when the 
solar activity in the southern hemisphere is dominant compared with when 
the solar activity in the northern hemisphere is dominant. Its reduction 
rates at four stations are all larger than those of the suppression due to other 
factors including the solar polarity effect on the GCR influx. We perform the 
student's t-test to see how significant these suppressions are. It is found that 
suppressions due to the solar North-South asymmetry as well as the solar 
polarity are significant and yet the suppressions associated with the former 
are larger and more significant. 
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12 1. Introduction 



It is well known that the level of galactic cosmic ray (GCR) influx to 
the terrestrial atmosphere is variable in time. Th e observed GCR influx is 



15 anti- correlated with the observed sunspot number 



Forbush 



1954 



Usoskin 



20081 ] . The higher level of the solar activity induces a stronger magnetic field 



strength in the heliosphere, and consequently reduces the amount of GCRs 
observed in the vicinity of the Earth. For instance, during the increasing 
phase of the solar cycle, the observed level of GCRs near the Earth has been 
gradually reduced. A long-term effect is also studied by Lockwood et al. 
(1999) using several geo- and solar-magnetic indices from 1964 to 1995 cor- 
responding to solar cycles from 20 to 22. They showed that measurements of 
the near-Earth interplanetary magnetic field reveal that the total magnetic 
flux leaving the Sun has risen by a factor of 1.4 since 1964, which interest- 
ingly reflects the increase of solar source magnetic flux, which also indirectly 
implies the decrease of GCR influx. 

In addition to the dependence of the GCR flux on the strength of the solar 
activity, the heliospheric c urrent sheet is asymmetric with respect to the heli 



29 ographic equatorial pla ne 



Munakata et al. 



Krymsky et al. 



2009 



Hiltula and Mursula 



2005 ; 



19831 ] . For instance, a southward offset of the current sheet 



3i has been reported while ULYSSES mon itored the cosmic-ray proton inten 



sity 



Smith et al. 



2000 



Simpson et al. 



rom September 1994 to May 1995 
19961 ] . We note that the observing epoch corresponds to the period when 

34 the solar southern hemisphere is magnetically dominant. As one of the solar 

35 northern- and southern-hemispheres is dominant in a given solar activity, 

36 the magnetic field of the heliosphere near the Sun is likely to follow a similar 



2 



37 behavior. IZhang et al. I |2002j have also reported that the electron temper- 

3B ature and magnetic field in the solar southern coronal hole are higher and 

39 stronger than those in the northern coronal hole during the period from 1992 

40 to 1997. This is also the time interval when the solar activity in the southern 

41 hemisphere is dominant, as we mentioned. 

42 In this paper, we investigate whether the level of the GCR influx shows 

43 any dependence associated with the solar North-South asymmetry using the 

44 neutron count rates obtained from 4 stations and sunspot data archived span- 

45 ning six solar cycles from 1953 to 2008. This is an interesting question not 

46 only in itself but also because of the following reason. According to what 

47 is mentioned above one may expect that the average location of the cur- 
4B rent sheet around the Earth can be quasi-periodically varying with the solar 

49 North-South asymmetry. Thus, the observed GCR influx may be changing in 

50 magnitude not only with the solar cycle but also with the solar North-South 

51 asymmetry. 



52 2. Data and Results 

53 The GCR data are recorded at ground-based neutron monitors, which 

54 detect variations in the low energy part of the primary cosmic ray spectrum. 

55 The lowest energy that can be detected at the top of the atmosphere depends 

56 on the geomagnetic latitude, and ranges from 0.01 GeV at stations near the 

57 geomagnetic poles to about 15 GeV near the geomagnetic equator. The 

58 National Geophysical Data Center(NGDC) tabulates the daily average GCR 
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59 count rate per hour for a number of stations. From the websitqj, we have 

60 taken the yearly mean of the GCR influx observed at 4 stations, where GCR 

61 data sets having a time-span longer than 20 years in both prior and post to ~ 

62 1980 are available: Moscow(55°N, 37°E, 2.42GV), Kiel(54°N, 10°E, 2.32GV), 

63 Climax(39°N, 106°W, 2.99GV) and Huancayo(12°S, 75°W, 12.92GV). 

64 We have adopted the yearly averaged wolf number (SN) from 1953 to 2008 

65 provided by Solar Influences Data Analysis Center (SIDcJ^I for the GCR-SN 
ee plot. We have also taken sunspot area data for the calculation of the solar 
67 North-South asymmetry from the NASA website! which is regularly updated 
es by Hathaway at the Marshall Space Flight Center (MSFC). Text files are there 

69 available containing the yearly averages of the daily sunspot areas (in units of 

70 millionths of a hemisphere) in the northern and solar southern hemispheres 

71 separately. 

72 In Figure 1 we show the running average of the solar North-South asym- 

73 metry, which is defined as the difference of the sunspot area between so- 

74 lar northern and southern hemispheres normalized by their sum, (An — 

75 As) /(Am + As), for the period of 1953 — 2008 as a function of time. In 

76 this plot the length of the moving window is 7 years. When values of the 

77 solar North-South asymmetry is positive the northern solar hemisphere is 

78 more active, and vice versa. Note that the dominance of the active hemi- 

79 sphere is reversed at ~ 1980. That is, in general the solar northern and 
so southern hemispheres are magnetically dominant before and after ~ 1980, 



1 http:/ /www. ngdc.noaa.gov/ 



http://sidc.oma.be/index.php3 



http : //solarscience. msfc.nasa.gov/greenwich.shtml 
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ai respectively. 

82 In Figure 2, we show yearly GCR fluxes at four neutron monitor stations 

83 as a function of the sunspot number. We divide the observed GCR flux and 

84 SN data sets into two sub-sets, respectively, separated at ~ 1980 for the 

85 present analysis. Open and filled circles represent the observed GCR count 
se rates in the sub-sets with time intervals before and after ~ 1980, respectively. 
87 Straight lines in each panel result from the least square fits of two sub-sets, 
ss It should be noted that the averaged value of the observed GCR count rates 

89 represented by filled circles is smaller than that represented by open circles 

90 while the slopes of the best fits in each panel are more or less the same. 

91 The mean differences of GCR count rates between two sub-sets in Moscow, 

92 Kiel, Climax and Huancayo are 2.83±0.33%, 2.50±0.44%, 2.72±0.18% and 

93 0.72±0.04%, respectively. What it means is that while the anticorrelation 

94 between the GCR flux and the sunspot number holds good regardless of ob- 

95 servation durations, the averaged GCR flux is apparently more reduced when 

96 the solar southern hemisphere is more active than the northern hemisphere. 

97 The relative difference shown is calculated from the best fit. We compute stu- 

98 dent's t and its probability for each station to see how significantly two mean 

99 values differ from. Student's t values for Moscow, Kiel, Climax and Huan- 

100 cayo are 9.50, 7,38, 6.28 and 5.21. Their false-alarm levels are all less than 

101 0.01. It should be noted that the percentage variations in absolute counting 

102 rates during the solar cycle are different for different geomagnetic latitudes. 

103 The effect is smaller at lower latitudes in agreement with the shielding effect 

104 of the Earth's magnetic field on high-energy charged particles. Note that 

105 the vertical rigidity cut-offs for Moscow, Kiel, Climax, and Huancayo are 
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106 2.42GV, 2.32GV, 2.99GV, and 12.92GV, respectively. As a result, the mean 

107 value of data points in each sub-set is considered as significantly different. 

108 This conclusion is found to be insensitive to the epoch when we cut the data 

109 set. We have repeated our analysis with a few more sub-sets separated at 
no different times rather than ~ 1980, and obtained similar results. 

in To compare our findings with results from other possible factors causing 

ii2 the suppression we have observed in the present work, we have carried on 

n3 our analysis with the data sets divided by three other criteria. By doing 

H4 so we may wish to demonstrate that the solar North-South asymmetry is a 

us more critical factor in reducing the observed GCR flux than any other solar 

lie variables compared in this paper. 

a? In Figure 3, we show the GCR influx observed at the Moscow station as 

us a function of the sunspot number, as an example. In panel (a), the data set 

n9 is divided by the level of the solar activity. That is, open and filled circles 

120 represent the observed GCR count rates in the period of strong and weak solar 

121 cycles, respectively. Open circles belong to the cycles 19, 21 and 22. Filled 

122 circles belong to the cycles 20 and 23. Note that cycles 19 and 20 correspond 

123 to the period prior to ~ 1980, and cycles 21, 22 and 23 to the period post to 

124 ~ 1980. The solar cycle is defined such that it begins and ends at the solar 

125 minimum. A clear difference in the mean value cannot be seen. In panel (b), 

126 we show the observed GCR influx as a function of the sunspot number with 

127 respect to the individual solar cycle. Different symbols represent different 

128 solar cycles as indicated by the number. Note that the solar cycle is defined 

129 such that it begins and ends at the solar maximum unlike the case in panel 

130 (a). Having picked up the cycles 20 and 23, and compared them with other 
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131 solar cycles by averaging, the results can be contrasted with panel (a). By 

132 defining the solar cycle in this way one may also see the polarity effect. In 

133 panel (c) the polarity effect on the observed GCR influx is shown. Open and 

134 filled circles represent the opposite polarity. T wo consecutive solar cycles 



135 have the opposite polarity. As reported earl ier 



1980 



Jokipii Ill989l ; ISabbah and Duldig 



)7] 



Nagashima and Morishita 



20071, we confirm that the amount 



137 of the observed GCR influx varies as the polarity of the solar magnetic field 

us is reversed. For comparison in magnitude, in panel (d) we show the result 

139 from the solar North-South asymmetry. It should be noted that the solar 

wo North-South asymmetry effect on the GCR influx is larger than any other 

mi effects. That is, the values of dl/I due to the solar North-South asymmetry 

142 and the polarity are 2.83±0.33 % and 1.36±0.50 %, respectively, as shown 

143 in Table 1. Similar conclusions can be drawn from other stations. 

144 It is also interesting, for instance, to note that the mean difference in the 

145 observed GCR influx is very large in cycles 20 and 23. These are the solar 
U6 cycles whose signs both in the solar polarity and in the solar asymmetry are 
147 simultaneously opposite. These two factors apparently compete with each 
us other so that they can be either amplified or compensated accordingly. Hence 
us we conclude that both effects should be considered with due care, particularly 

150 when one wishes to employ one of these effects into a sophisticated algorithm 

151 to predict anything accurately. 

152 In Table 1, to summarize results, we list the difference in mean values 

153 and results of the student's t-test for sub-sets divided according to the so- 

154 lar polarity and the solar North-South asymmetry. In the first and second 

155 columns, neutron monitor stations and the criteria in dividing the data sub- 
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156 sets are shown. In the third column the relative difference in the mean GCR 

157 flux is given. In the fourth, fifth, and sixth columns, results of student's 

158 t-test are given. We conclude that the suppression associated with the solar 

159 North-South asymmetry is larger and more significant than that with the 
wo solar polarity. We only consider in the present study the average long-term 
lei effect of the solar North-South asymmetry on the GCR influx. Though the 
162 solar North-South asymm etry also has sh ort er-term perio dicities than we re- 



163 ported in this paper[e.g., IChang 



2007a 



2008 



2009|| , those short-terms 



164 may not be notable due to the other effect such as a fluctuating shape of 

165 current sheet. 



we 3. Summary and Discussion 

167 We investigate whether the level of the GCR influx shows any dependence 

168 associated with the solar North-South asymmetry using neutron count rate 

169 data obtained from 4 stations and sunspot data archived spanning six solar 
no cycles from 1953 to 2008. We have found that the GCR influx observed 
ni at 4 neutron monitor stations are all suppressed when the solar southern 

172 hemisphere is more active. The amount of suppressions are 2.83±0.33% at 

173 Moscow, 2.50±0.44% at Kiel, 2.72±0.18% at Climax and 0.72±0.04% at 

174 Huancayo. These values are larger than the solar polarity effect observed 

175 at 4 stations (1.36±0.50% at Moscow, 1.67±0.57% at Kiel, 2.22±0.89% at 
ne Climax and 0.39±0.04% at Huancayo). We compute student's t values and 

177 its probability for each case, and come to conclusions that the solar North- 

178 South asymmetry effect on the GCR influx is larger than any other effects 

179 including the solar polarity effect. 
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Finally, we would like to add some speculation of its possible implica- 
tion on the terrestrial climate. As Svensmark and Friis-Christensen (1997) 
pointed out, the decrease in the observed GCR influx may cause the low 
terrestrial albedo through reducing low level cloud cover, and consequently 
lead to the higher heat content at the surface of the Earth. In a sense, 
therefore, the heat content at the surface of the Earth could be different 
from that deduced by the measured TSI alone [Goode et al. 1 120031 ] . We sus- 



pect that the net-radiative energy at the surface of the Earth should be 
determined by two parameters even in a simple climat e model based on an 



assumption that greenhouse fraction is constant [e.g., 



Dutton 



1995| . For 



instance, two parameters can b e the tota 
modulated with the solar cycle [Scafetta 



solar irradiance (TSI), wh ich is 



2009 and references therein!, and 



the Earth's albedo, which Svensmark and Friis-Christensen (1997) show vary- 
ing with the observed GCR influx modulated. Hence one may wish to re- 
examine the temperature anomaly as a function of solar activity including, 
for example, the solar North-South asymmetry since the correlation between 
the GCR i nflux and the sunspot number is subject to the sign of the sola r 



197 asymmetry 



Cho et al. 



2009 



Georgieva 



2002 



Georgieva et al. 



2005 



2007 ]. 



198 4. Acknowledgments 

199 We appreciate Katsuhide Marubashi for constructive comments on the 

200 GCRs distribution in the heliosphere and Kate Connors for careful reading 

201 the manuscript. We thank the anonymous referees for constructive comments 

202 which improve the original version of the manuscript. This work was sup- 

203 ported by the 'Development of Korean Space Weather Center' of KASI and 



9 



204 KASI basic research funds. HYC was supported by Basic Science Research 

205 Program through the National Research Foundation of Korea(NRF) funded 

206 by the Ministry of Education, Science and Technology (2009-0071263). 

207 References 

208 Chang, H.-Y., 2007a. Variation in North-South Asymmetry of Sunspot Area. 

209 Journal of Astronomy and Space Sciences 24 (2), 91-98 

210 Chang, H.-Y., 2007b. A New Method for North-South Asymmetry of Sunspot 

211 Area Analysis 24 (4), 261-268 

212 Chang, H.-Y., 2008. Stochastic properties in North-South asymmetry of 

213 sunspot area. New Astronomy 13 (4), 195-201 

214 Chang, H.-Y., 2009. Periodicity of North-South asymmetry of sunspot area 

215 revisited: Cepstrum analysis. New Astronomy 14 (2), 133-138 

216 Cho, I. H., Kwak, Y. S., Cho, K. S., Choi, H. S., and Chang, H.-Y., 2009. 

217 On the relation between the Sun and climate change with the solar North- 

218 South asymmetry. Journal of Astronomy and Space Sciences 26 (1), 25-30 

219 Dutton, J. A., 1995. Scientific priorities in global change research. In: Asrar 

220 and D. J. Dokken(Ed), The state of earth science from space. American 

221 Institute of Physics. 

222 Forbush, S. E., 1954. World-Wide Cosmic- Ray Variations, 1937-1952. Journal 

223 of Geophysical Research 59 (4), 525-542 



10 



224 Georgieva, K., 2002. Long-term changes in atmospheric circulation, earth 

225 rotation rate and North-South solar asymmetry. Physic and Chemistry of 

226 the Earth 27 (6-8), 433-440 

227 Georgieva, K., Kirov, B. and Bianchi, C, 2005. Long-term variations in 

228 the correlation between solar activity and climate. Memorie della Societa 

229 Astronomica Italiana 76 (4), 965 

230 Georgieva, K., Kirov, B. Tonev, P., Guineva, V., and Atanasov, D., 2007. 

231 Long-term variations in the correlation between NAO and solar activity: 

232 The importance of north south solar activity asymmetry for atmospheric 

233 circulation. Advances in Space Research 40 (7), 1152-1166 

234 Goode, P. R., Palle, E., Yurchyshyn, V., Qiu, J., Hickey, J., Rodriguez, P. 

235 Montanes, Chu, M.-C, Kolbe, E., Brown, C. T., and Koonin, S. E., 2003. 

236 Sunshine, earthshine and climate change: II. Solar origins of variations in 

237 the earth's albedo. Journal of the Korean Astronomical Sociery 36 (SI), 
S83-S91 

239 Hiltula, T. and Mursula, K., 2005. Heliospheric current sheet North-South 

240 asymmetry since 1926. In: D. Danesy, S. Poedts, A. De Groof and J. 

241 Andries(Ed), Proceedings of the 11th European Solar Physics Meeting 

242 " The Dynamic Sun: Challenges for Theory and Observations" , 600E, 123H 

243 Jokipii, J. R., 1989. The physics of cosmic-ray modulation. Advanced Space 

244 Research 9 (12), 105-1 19 

245 Krymsky, G. F., Krivoshapkin, P. A., Mamrukova, V. P., and Gerasimova, 



11 



246 S. K., 2009. North-south asymmetry of the heliosphere from cosmic-ray 

24? observations. Astronomy Letters 35 (5), 333-337 

248 Lockwood, M., Stamper, R., and Wild, M. N., 1999. A doubleing of the Sun's 

249 coronal magnetic field during the past 100 years. Nature 399, 437-439 

250 Munakata, Y., Hakamada, K., Mori, S., Yasue, S., and Ichinose, M., 1983. 

251 Cosmic ray North-South asymmetry related with the latitudinal angular 

252 distance of the earth from the heliospheric current sheet. Proceedings from 

253 the 18th International Cosmic Ray Conference 3, 358M 

254 Nagashima, K. and Morishita, I., 1980. Twenty- two year modulation of cos- 

255 mic rays associated with polarity reversal of polar magnetic field of the 

256 sun. Planetary and Space Science 28 (2), 195-205 

257 Sabbah, I. and Duldig, M. L., 2007. Solar polarity dependence of cosmic ray 

258 power spectra observed with Mawson underground muon telescopes. Solar 

259 Physics 243 (2), 231-235 

260 Scafetta, N., 2009. Total solar irradiance satellite composites and their phe- 

261 nomenological effect on climate, eprint Arxiv 0908.0792vl, submitted 

262 Simpson, J. A., Zhang, M., and Bame, S., 1996. A solar polar North-South 

263 asymmetry for cosmic-ray propagation on the heliosphere: The ULYSSES 

264 pole-to-pole rapid transit. Astrophysical Journal Letters 465, L69-L72 

265 Smith, E. J., Jokippi, J. R., Kota, J., Lepping, R. P., and Szabo, A., 2000. 

266 Evidence of a North-South asymmetry in the heliosphere associated with a 

267 southward displacement of the heliospheric current sheet. The Astrophys- 

268 ical Journal 533 (2), 1084-1089 

12 



269 Svensmark, H., Friis-Christensen, E., 1997. Variation of cosmic ray flux and 

270 global cloud coverage - a missing link in solar-climate relationships. Journal 

271 of Atmosphere and Solar- Terrestrial Physics 59 (11), 1225 - 1232 

272 Usoskin, I. G., 2008. A History of Solar Activity over Millennia. Living Re- 

273 views in Solar Physics 5 (3) 

274 Zhang, J., Woch, J., Solanki, S. K., and von Steiger, R., 2002. The sun 

275 at solar minimum: North-South asymmetry of the polar coronal holes. 

276 Geophysical Research Letters 29 (8), 77-1 



13 



1 .0 




-1.0 I , , I , , I 

1950 1960 1970 1980 1990 2000 2010 

Year 

Figure 1: The running average of the solar North- South asymmetry with a 7-year window. 
The vertical line is given to separate the data set into two sub-sets. Note that the solar 
northern- and southern- hemisphere is seen more active prior to ~ 1980 and post to ~ 
1980, respectively. Error bars are included. 
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Figure 2: Observed GCR count rates measured at Moscow(a), Kiel(b), Climax(c) and 
Huancayo(d) as a function of sunspot numbers. Open and filled circles represent the 
obseved GCR influx in the period before and after f980. Straight lines in each panel 
result from the least square fits of two sub-sets. 
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Figure 3: GCRs count rates measured at Moscow station as a function of sunspot numbers. 
Each straight line is the result of the least square fit. Different panels result from sub- 
sets divided by different criteria, i.e., (a) by the activity level of the solar cycles, (b) by 
the solar cycles, defined by the period between solar maxima, (c) by the solar magnetic 
polarity, (d) by the solar North-South asymmetry. 
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Table 1: Relative difference by the solar North-South asymmetry and polarity. 



Station Factors dl/I Student-t Deg. of freedom Prob.(%) 



Moscow 


Polarity 


1.36±0.50% 


2.83 


48 


>99.09 




Asymmetry 


2.83±0.33% 


9.50 


48 


>99.99 


Kiel 


Polarity 


1.67±0.57% 


3.89 


49 


>99.95 




Asymmetry 


2.50±0.44% 


7.38 


49 


>99.99 


Climax 


Polarity 


2.22±0.89% 


4.78 


51 


>99.99 




Asymmetry 


2.72±0.18% 


6.28 


51 


>99.99 


Huancayo 


Polarity 


0.39+0.04% 


2.21 


51 


>99.31 




Asymmetry 


0.72±0.04% 


5.21 


51 


>99.99 
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